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We first observe that the (co)domains of the g-deformed functions are some subsets of the 
(co)domains of their ordinary counterparts, thereby deeming the deformed functions to be incom- 
plete. In order to obtain a complete definition of g-generalized functions, we calculate the dual 
mapping function, which is found equal to the otherwise ad hoc duality relation between the ordi- 
nary and escort stationary distributions. Motivated by this fact, we show that the maximization of 
the Tsallis entropy with the complete g-logarithm and g-exponential implies the use of the ordinary 
probability distributions instead of escort distributions. Moreover, we demonstrate that even the 
escort stationary distributions can be obtained through the use of the ordinary averaging procedure 
if the argument of the g-exponential lies in (-oo, 0]. 
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^ ; I. INTRODUCTION 

Although Tsallis generalization of Boltzmann-Gibbs (BG) thermostatistics [l| has found many diverse fields of 
application there are some important issues in need of clarification such as the definition of temperature [1,0,[1,[1,0| 
and the associated central limit theorem 0, [jj [ToL ITTIl . One such subject of debate is the role of constraints in the 
maximization procedure of the Tsallis entropv~~[l2|. One choice of constraints consists of adopting the ordinary 
probability distributions, whereas a second choice (historically third choice) is adopted by considering the so-called 
escort distributions *2j] . Although many fundamental features of the underlying thermostatistics such as the Legendre 
structure seem to be preserved in both choices of constraints [12| . it is not easy to comprehend the need and the 
\ physical meaning of a new definition of probability distribution (and constraint) i.e., escort distribution. 
■ There have been many attempts to choose between these two distinct types of constraints. For example, one such 
0^ [ attempt has been made by considering the Shore- Johnson axioms in order to choose between these two probability 
distributions [llj]. However, it was later rejected on the ground of irrelevance Some recent attempts along these 
lines have been made concerning the stability of these distributions fl5j. It seems though that a final verdict is hard 
, to reach through these considerations, too [la [TtJ • 

However, all these aforementioned attempts explicitly depend on the g-deformed functions present in the definition 
of the Tsallis entropy expression. In fact, the mathematical structure of the g-thermostatistics is usually obtained by 
replacing the ordinary logarithmic and exponential definitions by their g-deformed counterparts. On the other hand, 
these deformed generalizations are not complete as their ordinary counterparts i.e., the deformed (co)domains do not 
extend to the whole appropriate range and are dependent on the deformation parameter q. This physically implies that 
these deformed functions can map the arguments only to a restricted codomain, not due to the constraints imposed by 
the system, but due to the incompleteness of the deformed functions. Moreover, due to the feature of incompleteness, 
these deformed functions are not invertible in entire (co)domain. In other words, the g-thermostatistics based on the 
incomplete deformed functions is unable to explain the asymptotic inverse power law decays in a consistent manner. 
Therefore, all discussions based on these incomplete deformed functions are in difficulty of being incomplete, too. This 
is also the case when one tries to discuss the nature of constraints in the maximization procedure of the associated 
Tsallis entropy. 

In order to overcome the incompleteness of the g-deformed functions, we propose the existence of a dual mapping 
d(q) in such a way that the deformed functions become invertible in the entire domains of their ordinary counterparts. 
We also note that the criterion of completeness based on finding a dual function d(q) is general (see Ref. [18| for 
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details) and can be applied to other deformed functions such as the ones by Kaniadakis [l9[ and Abe [20| . Doing this, 
we also eliminate the dependence of the intervals on the deformation parameter q. This enables the only constraints 
to be the ones imposed by the system. Then, we maximize the Tsallis entropy by using the complete g-deformed 
functions. The main result of this complete maximization procedure is to show that the escort distributions used so 
far are merely a means to compensate for the incompleteness of the underlying mathematical structure, and hence 
irrelevant to the formalism, when one uses the complete deformed functions. The ordinary probability distribution 
does not only suffice to maximize the Tsallis entropy, but also proves to be the correct constraint to be used in the 
maximization procedure. We also show how the escort distribution emerges through ordinary averaging procedure, 
and comment on its relevance in a consistent framework. 

The paper is organized as follows. In Section II, we first demonstrate the incompleteness of the g-deformcd functions 
and then define their complete counterparts. We maximize the Tsallis entropy based on complete deformed functions 
to show that the escort distributions are redundant in Section III. Concluding remarks is presented in Section IV. 



II. COMPLETE g-GENERALIZED FUNCTIONS 



Tsallis generalization is based on the following definition of deformed functions 

x 1 - 11 - 1 



/»(*) == 



h q {x) := l + (l-q)x 



(1) 



with h q — f„ , fq—yq =x(x) — ln(x) and h q —, qo =i(x) — exp(ir). The aforementioned functions are characterized by the 
following identities 



hq{-x)hq(x) 

I 1/(1-9) 



A, (a:), 



(2) 



with e q (x) := x "+* q q - 2 (x > 0) and A, (a;) := [l - (1 - qfx 2 ]^ 1 q) where [X}+ = max{A,0}. The functions 
Q q and A 9 are related in the sense that ® q (x) ^ <^> A q (x) ^ 1 or equivalently Q q (x) = A q (x) = 1. 

Before proceeding further, we define the following intervals Lo := {x G R : < x ^ 1}, Li := {x G R : x 5* 1} 
so that LoljLi = R + := {x G R : x > 0} i.e., the set of all positive real numbers. We also have the intervals 
R„ := {x e R : x ^ 0} and R^ := {x S R : x ^ 0}, respectively so that R = Rq [j I 
Concerning the domain and codomain of the g-functions in Eq. ([!]) , we observe 



Iq is the set of all real numbers. 
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Eq. ([3]) shows that the functions f q and h q represent a possible generalization only in some subsets of the codomains 
of the ordinary logarithmic and exponential functions when q ^ 1. In other words, these g-deformed functions are not 
invertible in the entire (co)domains of their ordinary counterparts. Since they are not generally invertible, a genuine 
q-exponential decay is ill-defined. Furthermore, from Eq. ©, we see that the function h q (—x) does not describe 
a g-exponential decay l/h q (x), as has been frequently considered in literature. Therefore, we call these functions 
as incomplete generalized functions. Searching for the origin of the aforementioned incompleteness we inspect e.g., 
the ordinary logarithm (in natural base e). We observe that In : R + — > R has an argument intercept at the point 
x = 1, dividing the logarithmic domain R + into two subdomains Lo and Li. The relation between the images of the 
log-arguments x £ Li and 1/x G L is represented by 



f qo (x)+f qo (l/x) = 



(4) 



Since the right hand side of Eq. (|4|) is equal to zero, the logarithmic codomain for x G R + is R. If the sum of the 
above logarithmic terms would be equal to Q(x), <d being a real non-singular function, then the codomain of the 
logarithmic function would depend on the codomain of Q(x) as well i.e., In : R + -tMCI, where M is some subset 
of the set of all real numbers. Moreover, inverting Eq. one would obtain exp(— x) ^ l/exp(a;) for Q(x) ^ 0. It 
becomes then evident that the origin of the incompleteness lies on the existence of the functions Q q and K q . Due to 
this incompleteness, g-thermostatistics is bound to accept arguments and yield results only in limited (co)domains. 



Furthermore, in these (co)domains the deformed exponential function 
negative arguments. 



h q exhibits different features for positive or 
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In order to define complete g-generalized functions, one must try to equate ® q (x) to zero (or equivalent A q to 
unity). Before proceeding further, let us denote the range of validity for the parameter q by A q C R. A boundary 
value of A q is always qo = 1 e.g., A q := (a, 1] with a < 1 or A q := [1, a) with a > 1. For defining complete deformed 
functions, we assume the existence of a function d(q) : A q — > B q (B q C R) such that 

f q (x) + f d{q) (l/x) = f q (l) (5) 

with lim g ^i d(q) = 1. Setting x = 1 in Eq. ([5]), we obtain = in accordance with the ordinary logarithm. 

It can be seen that the above relation is valid for any set of parameters by substituting q' = d(q) into the above 
equation. Therefore, Eq. ^ can be rewritten as 

f q (x)=-f d(q) (l/x) (6) 

in analogy to Eq. (|4]). The mapping d(q) is called as the dual function and the correspondence q <-> d(q) as duality 
relation. 

We further observe that the argument x = 1, since = 0, divides the (co)domain R + (R) into the following 

sub (co) domains R + = Lo IjLi (R = R~ U^cf )■ Consequently, we are led to the complete q-deformed definitions 

!fd( q ) ■ Lo — > R I h d ( q ) : R — > L 

and exp ? = In" 1 := < , (7) 

f q :h 1 ^R+ [h q :R+^U 

with q G A q . We can then, for the complete generalized logarithm and exponential given above, verify that 

ln q (l/x) = fd{ q ){l/x) = -f q (x) = -ln 9 0), 

for x £ Li and x £ Rj, respectively. 

The parameter range A q can be determined by requiring the fulfillment of the following limits satisfied by the 
ordinary functions 

lim cxp (x) = 1, lim \n q {x) = 0, (8a) 



lim exp (x) = 0, lim ln„(ir) = — oo, (8b) 

x — > — oo y x — >0 



x^l 

lim exp Ax) — oo, lim ln„(a;) = oo, (8c) 

X — >OG ^ X — >OG 

lim -^-exp (x) = 0, lim ln„(rc) = 0. (8d) 

x^ — oo dx x^oo dx 

Condition (|5aj) ensures the continuity of ln q (x) and exp q (x) at the points x = 1 and x = 0, respectively. The 
behavior of the complete deformed functions at the boundaries of their domains is determined from Eqs. (|8bp and 
(I8cl) . The condition (|8d[) allows one to preserve the same absolute maximum and minimum values of the ordinary 
functions. 

Having obtained the criterion of completeness in detail, we see that the first step is to calculate the dual function 
d(q). We can calculate it from Eqs. |T]) and © as 

d(q) = 2-q. (9) 
The requirements listed in Eq. ([8]) confines the values of the deformation parameter q into the following interval 

A- (o,i]. (io) 

The image of A q under the dual mapping d(q) i.e., B q can be calculated from Eqs. ([9]) and (fT0|) as 



B q := [1,2). 



(11) 
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Having explicitly obtained the dual mapping function d(q) and the range of validity of the deformation parameter 
q G A q , we can now write the analytical expression of the complete ^-generalized functions 



r xi- 1 - i 



big (a:) := < 



3-1 

x x -i - 1 



X £ hi 



exp Ax) 




, x eR 



, a; <E 



(12) 



in accordance with Eq. ([7]). We note that a slightly different expression for exp q is also obtained in Ref. [2l| in the 
context of cut-off prescriptions associated with the ^-generalized exponential. The division of the (co)domains yielding 
to a complete definition of the deformed functions was first noticed in Ref. [21]. However, this observation has not 
been pursued further as a criterion of completeness therein. Furthermore, the current definition is in accordance with 
the trace- form entropic definition based on f q (see next Section). 

The summary of the above results can be provided in a rather compact form below, considering two parameter 
ranges, A q and B q , and keeping a single expression in the entire (co)domains instead i.e., 



i_ — f . x£L I I seLi 
m 9 — Jq ■ qe[l,2)U 96(0,1] 



96] 



l,2)UoG(0,°ll' eXP 9 



xeR I I 
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eL I I zeL 
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(13) 



III. THE MAXIMIZATION OF THE TSALLIS ENTROPY: ORDINARY VERSUS ESCORT 

DISTRIBUTIONS 

The Tsallis entropy expression reads 



Sq(Pi) 



(14) 



where Pi is the probability of the ith-configuration and Q q is the maximum configuration function of the respective 
g-ensemble. The above expression can be rewritten in terms of the deformed g-logarithm given by Eq. (JTJ) as 



S q(?i) = ^Pifq (1M) • 



(15) 



In the literature of g-thermostatistics, there is no consensus on how the Tsallis entropy must be maximized. There 
are two distinct choices for the constraints used throughout the literature. The first one is based on considering the 
ordinary probability distributions when one averages the constraints in the functional to be maximized. This can be 
written as 



S^S q {pi)-aJ2Pi-Pj2P^ =0 - 



The maximization of this functional yields the stationary distributions, denoted by pi, of the form 



(16) 



Pi 



h q {ei) 



(17) 



Second choice rests on the definition of escort distributions. The escort distributions P, are defined as 



k 

The maximization of the Tsallis entropy with the escort distributions becomes 



p '-m (18) 
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(19) 



The above maximization yields the following stationary distribution 



- t 7~: = h q (-ei). (20) 

h2-q(£i) 

However, both of these results rely on the use of Tsallis entropy S q and this entropy in turn is based on the 
incomplete g-logarithm. Considering that one uses the complete deformed functions described above, we can have a 
more consistent look into the nature of stationary distributions obtained in the maximization procedure. 

The first observation considers the argument of the (/-logarithm used in the definition of the Tsallis entropy. The 
microstate probabilities vary between and 1 i.e., p, e [0, 1] so that 1/pt will be equal to or greater than 1. This 
means that the argument of the complete q- logarithm takes values in Lj . Then, the associated generalized statistics 
is characterized by the deformation parameter range A q = (0,1] according to Eq. (| 13|) . This parameter range 
corresponds to the g-exponential from Rq to Li , implying that the respective generalized exponential decay is described 
by the function l/h q (x). Inspecting Eq. (fT7|) , we see that this corresponds to the stationary distribution obtained 
from the maximization of the Tsallis entropy with ordinary constraints. 

Since the above considerations show that the maximization of Tsallis entropy must be carried out with ordinary 
probability definitions, one might wonder how the escort distributions emerge at all. In order to shed light on this 
issue, we inspect Eqs. (TT21) . (fl3l) and (|2"0|) to see that the stationary probability distribution obtained from the 
constraints averaged with the escort distribution corresponds to the parameter validity range q £ [1,2) and thus it 
should be related to the expressions transformed under the mapping d(q) E B q . This parameter range corresponds 
to the complete (/-exponential when x e i.e., when the argument of the (/-exponential equal to or less than 0. 
Moreover, if we write the Tsallis entropy compatible with this range i.e., 

s d{q) {Pi) = £ Pih {q) {i/Pi) = - ]T p*fM = E Pi a _ 1 ( 21 ) 

i=l i=l i=l * 

and normalize it with ordinary constraints, we see that it yields the stationary distribution given by Eq. I|20p i.e., 
Pi = I//12— q( s i)- This is the stationary solution one obtains using the constraints averaged with escort distributions. 

To summarize, both types of stationary distributions, ordinary and escort, can obtained through the use of ordinary 
constraints. However, it is the stationary distribution obtained through the ordinary averages, which is consistent 
with the range of the arguments, since the argument of the (/-logarithm i.e., 1/pi, is equal to or greater than 1. In 
other words, the escort distributions are only ad hoc means of conforming the incomplete (/-thermostatistics to the 
parameter range lying between 1 and 2. Moreover, it is worth mentioning that this result is in agreement with the 
concept of the entropy extensivity. Indeed, considering Tsallis entropy in Eq. (|15[) for equal probabilities pi — l/Q q 
and demanding the fulfillment of the extensivity property with respect to the variable £i we are able to determine the 
maximum configuration function given by Cl q (£i) = h q (ei) => pi = V^g ( £ i)- 

In contrast to the incomplete definitions, both ordinary and escort stationary distributions becomes identical in 
the context of the complete definitions presented herein, since exp 9 (— x) — 1/ exp ? (x), and can be obtained through 
ordinary averaging procedure. 



IV. CONCLUSIONS 



The g-thermostatistics is fundamentally based on q-deformed functions i.e., g-logarithm f q and ^-exponential h q . 
Therefore, it is essential to use valid definitions of these generalized functions. However, the aforementioned func- 
tions are not bijective in the entire (co)domain of their ordinary counterparts, implying that they do not represent 
complete generalizations, since they are not be invertible in the entire (co)domain. This incompleteness creates math- 
ematical (and physical) discrepancies, which makes the use of primary g-functions deficient. On the other hand, the 
maximization procedure of Tsallis entropy necessitates the use of the (/-logarithm in its definition and its inverse 
i.e., (/-exponential in obtaining the stationary solution resulting from this maximization procedure. Therefore, it is 
important to remaximize the Tsallis entropy with the corresponding complete definitions of the generalized functions. 
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The maximization of Tsallis entropy in terms of the complete g-generalized functions shows that the correct station- 
ary solution is the one associated with the ordinary constraint instead of the widely used escort stationary distributions. 
In this sense, the definition of the escort distributions is redundant. 

It is worth noting that the use of (/-generalized functions f q and h q is not solely limited to the entropy maximization 
procedure. The inspection of these functions in a more "complete" framework shows that the escort distributions 
are nothing but ad hoc means of extending the so far incomplete definitions of the deformed functions to the whole 
range of the deformation parameter. In other words, the escort distributions become present whenever the incomplete 
primary deformed functions cannot map certain regions in the deformation parameter space q, namely, q S [1,2). 
This region is physically important. Nevertheless, once we define the complete g-deformed functions, all ranges of 
the deformation parameter can be accessed in an invertible manner by these complete functions, thereby showing the 
redundancy of the escort distributions in g-thermostatistics. 

Last but not least, all the definitions of the q-thermostatistics based on the deformed functions must be revised 
accordingly in terms of complete deformed functions. One such important example is the recently defined q-Fourier 
transform 22], since it is closely related to obtaining a (possible) central limit theorem, whose basin of attraction 
is g-Gaussian. We finally note that the recently developed maximization procedure based on the definition of the 
q- Fourier transform must be revised accordingly too in order to assess the true role of the escort distributions (23j . 



Acknowledgments 

We thank U. Tirnakli for a careful reading of the manuscript. GBB was supported by TUBITAK (Turkish Agency) 
under the Research Project number 108T013. 



[1] C. Tsallis, Possible Generalization of Boltzmann-Gibbs Statistics, J. Stat. Phys. 521/2 (1988) 479. 

[2] C. Tsallis, Introduction to onextensive statistical mechanics: approaching a complex world, New York: Springer, 2009. 

[3] S. Abe, Temperature of nonextensive systems: Tsallis entropy as Clausius entropy, Physica A 368 (2006) 430. 

[4] S. Abe, Heat and entropy in nonextensive thermodynamics: transmutation from Tsallis theory to Renyi- entropy-based 

theory, Physica A 300 (2001) 417. 
[5] S. Abe, General pseudoadditivity of composable entropy prescribed by the existence of equilibrium, Phys. Rev. E 63 (2001) 

061105. 

[6] S. Abe, S. Martinez, F. Pennini and A. Plastino, Nonextensive thermodynamic relations, Phys. Lett. A 281 (2001) 126. 
[7] H. Hasegawa, Nonextensive thermodynamics of the two-site Hubbard model, Physica A 351 (2005) 273. 
[8] U. Tirnakli, C. Beck and C. Tsallis, Central limit behavior of deterministic dynamical systems, Phys. Rev. E 75 (2007) 
040106. 

[9] U. Tirnakli, C. Tsallis and C. Beck, A closer look at time averages of the logistic map at the edge of chaos, Phys. Rev. E 
79 (2009) 056209. 

[10] P. Grassberger, Proposed central limit behavior in deterministic dynamical systems, Phys. Rev. E 79, 057201 (2009). 
[11] A. Pluchino, A. Rapisarda and C. Tsallis, Nonergodicity and central limit behavior in long-range Hamiltonians, EPL 80, 
26002 (2007). 

[12] C. Tsallis, R. S. Mendes, and A. R. Plastino, The role of constraints within generalized nonextensive statistics, Physica A 
261, (1998) 534. 

[13] S. Abe, G. B. Bagci, Necessity of q- expectation value in nonextensive statistical mechanics, Phys. Rev. E 71, 016139 (2005). 
[14] S. Abe, Generalized molecular chaos hypothesis and H-theorem: Problem of constraints and amendment of nonextensive 

statistical mechanics, Phys. Rev. E 79, 041116 (2009). 
[15] R. Hanel, S. Thurner and C. Tsallis, On the robustness of q- expectation values and Renyi entropy, EPL 85, (2009) 20005. 
[16] S. Abe, Anomalous behavior of q-averages in nonextensive statistical mechanics, accepted to J. Stat. Mech. (2009), 

larXiv:0906.2908l 

[17] J. F. Lutsko, J. P. Boon and P. Grosfils, Is the Tsallis entropy stable?, EPL 86, (2009) 40005. 

[18] Th. Oikonomou and G. B. Bagci, Complete versus incomplete definitions of the deformed logarithmic and exponential 
function submitted (2009). 

[19] G. Kaniadakis, M. Lissia and A.M. Scarfone, Two-parameter deformations of logarithm, exponential, and entropy: A 
consistent framework for generalized statistical mechanics, Phys. Rev. E 71, 046128 (2005). 

[20] S. Abe, A note on the q- deformation-theoretic aspect of the generalized entropies in nonextensive physics, Phys. Lett. A 
224 (1997) 326. 

[21] A. M. Teweldeberhan, A. R. Plastino and H. G. Miller, On the cut-off prescriptions associated with power-law generalized 

thermostatistics, Phys. Lett. A 343 (2005) 71. 
[22] C. Tsallis, A. R. Plastino, and R. F. Alvarez-Estrada, Escort mean values and the characterization of power-law-decaying 

probability densities, J. Math. Phys. 50, (2009) 043303. 



7 



[23] G. B. Bagci and U. Tirnakli, On the way towards a generalized entropy maximization procedure, Phys. Lett. A, in press 
(2009). 



